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Abstract. - The noncentrosymmetric CIO4 anions of the organic superconductor (TMTSF)2C104 
order below 24K. The size of domains where the anions are ordered is substantially dependent 
on the cooling rate which is a key parameter for the stability of the low temperature electronic 
ground states. We study the effect of the cooling rate on the SC phase within a self consistent 
approach in the framework of the time dependent Ginzburg-Landau theory taking into account 
the superconducting fluctuations. We derive the superconducting transition temperature which is 
found to decrease with increasing cooling rate in agreement with recent experimental data. 



The Bechgaard salts (TMTSF) 2 X, where TMTSF de- 
notes the organic molecule tetramethyltetraselenafulva- 
lene and X stands for the inorganic anion (X= VFq, CIO4, 
Re04 etc.), continue to be under close scrutiny due to their 
rich phase diagram exhibiting a variety of ground states 
such as metallic and superconducting (SC) states as well 
as spin density wave (SDW) phases [1,2]. A key feature 
of this phase diagram is the proximity of the insulating 
SDW phase to the superconducting one, which raises up a 
worthwile question about the interplay between these two 
states: wether they coexist or compete. 
Transport measurements on (TMTSF) 2 PF 6 [3,4] revealed 
the presence of an inhomogeneous superconducting state 
in the vicinity of the border between the SDW phase 
and the superconducting one. This feature has been as- 
cribed to the coexistence of macroscopic SDW and super- 
conducting slabs under high pressure. Contraray to the 
(TMTSF) 2 PF 6 salt, (TMTSF) 2 C10 4 undergoes a super- 
conducting phase transition at ambient pressure. The sta- 
bility of this phase is substantially dependent on the rate 
at which the sample is cooled through the structural tran- 
sition temperature Tao Below this temperature, the non- 
centrosymmetric CIO4 anions order by alternating their 
orientations along the b direction perpendicular to the 
stacks of the TMTSF molecules arranged along the a axis. 
In the relaxed state, obtained by slowly cooling the sam- 
ple at a rate less than O.lK/mn, the CIO4 anions ordering 
occurs at 24K at ambient pressure and superconductivity 
appears at 1.2K. However, in the quenched sample, rapidly 



cooled at a rate of at least 50K/mn, the anion ordering can 
no more be achieved and the anions are randomly orien- 
tated. In this case, the superconducting ground state is 
destroyed and a SDW phase is established at 6K. 

What is remarkable with the anion ordering (AO) is the 
possibility to have intermediate states with a mixture of 
SDW and SC regions. This feature has been first reported 
by transport [5] and specific heat [6] measurements at dif- 
ferent cooling rates which have been corroborated by X 
ray diffraction studies [7] . Based on a muon spin rotation 
study, Greer et al. have also observed, for intermediate 
cooling rates, a mixed state with macroscopic sized do- 
mains of superconducting and SDW phases [8]. 

Moreover, it has been found that, by increasing the cool- 
ing rate, the SC transition temperature decreases and ends 
at 0.7K for a cooling rate of 36K/min according to Meiss- 
ner signal measurements [9,10]. 

The cooling rate dependence of the SC transition tem- 
perature T c has been recently re-investigated by Joo et al. 
who derived the phase diagram of (TMTSF) 2 C104 for dif- 
ferent cooling rates [11,12]. The results revealed that T c 
starts to decrease rapidly for small quenching rates (less 
than 12K/mn) but the decrease is slowed down as the cool- 
ing rate increases and T c tends to saturate at 0.9K before 
the suppression of the SC phase at a rate of 14K/mn. 

For the rates from HK/mn to 14K/mn, the authors re- 
ported the coexistence of a SC state and a SDW phase. 
The latter is revealed by the upturn of the resistivity pre- 
ceding its drop at the SC transition temperature in accor- 
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dance with the early measurements of Schwenk et al. [5]. 

This behavior has been ascribed to the presence of or- 
dered CIO4 anions domains where superconductivity de- 
velops and disordered regions where SDW phase appears. 
These two ground states do not coexist microscopically 
but are rather segregated as pointed out in Ref. [5] . 

For slow cooling rates, the size of disordered domains is 
small in such a way that they can be considered as defects 
or non-magnetic scattering centers, which are known to 
induce a sharp decrease of the transition temperature of 
non-s wave SC phase according to the Abrikosov-Gor'kov 
theory [13]. Therefore, the sensitivity of T c to the cooling 
rate in the (TMTSF)2C1C>4 has been considered as a proof 
for an unconventional SC state [10,12]. 

As the cooling rate increases, the disordered domains 
can no more act as point defects and the sensitivity of T c 
to the quenching is then more and more reduced which ex- 
plains the slowing down of the rate at which T c decreases 
as found by Joo et al. [11, 12]. 

In the quenched state, the SDW islands merge to form 
larger domains at the expense of the ordered ones. As a 
result the SC phase is completely suppressed and a pure 
SDW phase is stabilized. It is worth to note that a pure 
SDW (SC phase) has not been observed, there are always 
tiny ordered (disordered) domains where the SC (SDW) 
phase develops. 



A first attempts to interpret the effect of the AO on the 
ground state of the (TMTSF) 2 C10 4 salt was undertook 
by Hcriticr et al. [14] who derived a rough estimate of the 
SDW transition temperature in the relaxed state and in 
the quenched one. The authors argued that AO may per- 
turb the nesting properties of the Fermi surface and as a 
consequence, the SDW phase is destabilized at ambient 
pressure in the relaxed state. However, the possibility of a 
phase mixture has not been addressed. In Ref [15], renor- 
malization group calculations were proposed taking into 
account the anion ordering which is assumed to develop 
over the whole sample. The results showed that the sin- 
glet SC state is the most stable phase in the relaxed state 
while the SDW phase will develop in the quenched state. 
For the intermediate cooling rates, the results revealed a 
microscopically coexistence between the SC and the SDW 
phases while experiments suggest that these phases are 
rather segregated as discussed above. To deal with this 
phase segregation, a somewhat different renormalization 
group (RG) approach has been proposed based on a struc- 
ture of ordered and disordered anion domains. According 
to this model, the segregated phase is expected to be the 
unique possible description of the intermediate state [16]. 
However, these RG calculations are carried out in the one 
dimensional regime far from the low temperature ordered 
states. 

As far as we know, no theoretical study of these states, 
taking into account AO, has yet been provided, in partic- 
ular for the intermediate cooling rate. 



In this letter, we propose a model to describe the 
dependence of the SC transition temperature on the 
cooling rate within a selfconsistent approach based on 
the time dependent Ginzburg Landau (TDGL) theory 
and taking into account the superconducting fluctuations 
which cannot be neglected in Bcchgaard salts as argued 
experimentally [11,12]. This issue, which is the keystone 
of the present work, has not been addressed so far. The 
results deduced from our model are consistent with recent 
experimental data. 

We assume that the (TMTSF)2C104 sample is formed 
by ordered domains with a SC phase, embedded in 
a matrix of insulating islands where CIO4 anions are 
disordered. We take for simplicity a layered structure 
where in each layer, corresponding to the most conducting 
(a, b) plane, we consider an array of identical ordered 
domains separated by disordered ones which are also 
assumed to be identical. The tunneling between the 
ordered superconducting domains is insured by Josephson 
couplings denoted Ji and J2 in the a and the b direction 
respectively while the interplane Josephson coupling is 
parametrized by Jo- 

The average size of the SC island is assumed 
to be greater than the in-plane coherence length 
Zo = U = vCT-sooA [7]. 

To bring out the importance of the SC fluctuations 
in (TMTSF)2C104, we evaluate the Ginzburg number 
Gi which is the fundamental parameter governing the 
strength of the fluctuations. G;, indicating the tempera- 
ture range around the critical temperature T c where fluc- 
tuations are not negligible, is given by: 



Gi = 



ST 
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here 7 = is the anisotropy parameter where £ a b and 
£ c are the coherence lengths in the (a, b) plane and along 
the c axis respectively. 0o is the quantum flux and A is 
the penetration depth. 

Assuming £ a b =~ 500A, 7 ~ 20 and A = 17/xm [1], 
we found Gi ~ 0.1, which is considerably large as in 
high-T c . However, in conventional superconductors the 
Ginzburg number is small (Gi ~ 10~ 8 ). The large value 
of G l in (TMTSF) 2 C10 4 means that the SC fluctuations 
spread over a wide interval of temperature compared to 
the critical temperature T c , which is consistent with the 
departure from a Fermi liquid behavior of the temperature 
dependence of the c axis resistivity measured around 5 K 
and down to T c ~ 1.2 K [11, 12]. 

Due to the superconducting fluctuations in the granular 
structure considered in our model, we expect deviations 
from the mean field theory value To corresponding to the 
transition from the metallic state to the pure supercon- 
ducting phase where all anions are ordered. We will take 



Tn = T, 



exp 



Gi T, 



exp 



1.8 K where T, 



exp 



1.5 K is the 
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experimental SC transition temperature found recently in 
a very slowly cooled pure sample [17, 18]. 

To study the stability of the SC phase, we compare its 
Ginzburg-Landau free energy F s to that of the normal 
state F norrn [19]. 
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where X (Y) is the intra-domain coordinate along the a 
(b) axis while l\ and 1% are respectively the lengths of the 
SC domain along the a and the b directions. ip n ij is the SC 
order parameter in the n th plane and in the domain 
where i (j) labels the SC island in the a (b) direction. The 
a coefficient is given by: a — a e where ao = 7i 2 /2to£q 
and e = ln(T/T ) while b = ^ K 2 e^ 2 /2m 2 , k is the GL 
parameter, eo = 2e is the pair electric charge and to is the 
effective pair mass in the (ab) plane. 
The coupling parameters are written as: 



Ji = 



2m*l'{ 



J 2 = 



2m* l'i 



and Jq 



2m c s 



(2) 

with to* and m c are the effective pair masses in the SC 
domain and along the c axis. We take to* = m and we 
introduce the anisotropic parameter 7 = £0 /£ c = y/m c /m, 
where £ c is the coherence length is the c direction. l[ and 
V 2 are respectively the size of the insulating domains in the 
a and b directions whereas s is the intcrplanc distance 

The dynamics of the SC order parameter satisfies the 
time dependent Ginzburg-Landau (TDGL) equation: 



at 



OF 
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Here T^ 1 = ttH 3 /16to^q/cbT is the relaxation rate of the 
order parameter whereas C, n ij {t, t) are the Langevin forces 
governing the thermodynamical fluctuations and which 
fulfill the Gaussian white-noise law [19]: 

(U'(r,i)fe(r '.*')) = 2T^k B TS(f- f')S(t t') 5 -^f- 

with f = (X + i(h + l[),Y + j(l 2 + l' 2 ),ns) and 
r' = X + i'(h + Y + j'(l 2 + l' 2 ),n's). 

Taking the derivative of the free energy F (Eq.l) with 
respect to tp nij and replacing b\ip n ij\ 2 ip ni j by b(\ip ni j\ 2 }ip ni j 
to get a linear problem, the TDGL equation writes as: 



dt 



+ b(\lpnij\ 2 )^nij + J\ (2lpnij - V'ni+l j - lpni-1 j) 
+ J2 {2lpnij — 1pnij+\ ~ i>nij-l) 

+ Jo (Zlpnij ~ ipn+lij - fpn-lij) (3) 



We normalize the reduced temperature e: 



= e+~(\<Pnij\ 2 ) 



(4) 



and {\ip n 



will be determined self-consistently. We 



also define the Fourier transform of ip n ij as 

d 2 k f~ dq 



where 
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Equation 3 takes then the form: 

Cqi^xj kyi — 



„_! d h 2 kl n 2 k 2 v 



2m 2m 

+ 2J 1 (l-cos(k x (h + l' 1 )) 
+ 2J 2 (1 -cos(k y {l 2 +1' 2 )) 
+ 2J (1 - cos(qs)]ip q (k x ,k y ,q,t) (5) 

with 

(C q (k,t)Q,(k',t')) =2 r^k b T(27T) 3 S(k-k / )5(q-q') 
x 6(t-t') 

where k = (k x , k y ) and a = a + b(\ip n ij\ 2 ). 

Using the Green function method [19], one can solve 
Eq.5 and determine {\ipnij\ 2 ) by straightforward calcula- 
tions [20]. 

The self-consistent equation Eq.4 can now be solved to 
derive the critical temperature T c at which e = 0. 

In figure 1 we have depicted the superconducting transi- 
tion temperature as a function of the inverse of the reduced 
coupling J[ = J1/00 along the a axis, which is expected 
to decrease by quenching the sample since the distances, 
l[ and 1' 2 , separating two SC domains increase with cool- 
ing rate, putting at disadvantage the Joscphson tunnel- 
ing. Therefore the 1/ J{ parameter mimics the cooling 
rate. The exact dependence of the Josephson coupling on 
the cooling rate may be deduced from the behavior of the 
inter-domain distances as given by Eq.2. However, such 
behavior has not been determined experimentally since X 
ray diffraction measurements reported only the size pa- 
rameters of the ordered domains (li and I2) which are 
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found to be reduced by quenching [7]. We, then, assume 
that l[ and 1' 2 increase with the cooling rate and as a con- 
sequence 1/ Jj will be enhanced by quenching. 



1.5 - 
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Fig. 1: Superconducting transition temperature of the 
(TMTSF)2C104 as a function of the inverse of the reduced 
Josephson coupling along the a axis which mimics the cooling 
rate. Numerical results are obtained for h = li = 800A, where 
h and h are the lengths of the SC domain along the a and the 
b directions respectively. The arrow indicates the mean field 
theory value of the SC transition temperature To. 

As shown in Fig.l, T c , which is lower than the mean field 
theory value T due the SC fluctuations, decreases rapidly 
for small cooling rates for which the SDW domains can be 
regarded as point defects having drastic effect on a non- 
s wave SC phase. For larger cooling rates, T c decreases 
smoothly and tends to saturate with further quenching 
before ending at a critical value of the Josephson coupling 
J[ [10, 11]. It is worth to stress that, actually, T c does not 
saturate and continue to decrease but with a very small 
rate. However, the greater the cooling rate, the larger 
the SDW domains, the smaller the anion gap opening on 
the Fermi surface. The system, which is characterized by 
a four sheet Fermi surface in the relaxed state, is then 
more and more close to a one band Fermi surface as in 
(TMTSF) 2 PF 6 . Therefore, the nesting properties are im- 
proved and the transition temperature from the metallic 
state the SDW phase will then increase by quenching [14]. 

We have looked for the effect of the SC domain size on 
T c since the domains are expected to shrink by quench- 
ing [7]. However, we found that T c is almost independent 
of the size parameters l\ and h of the ordered domains. 
This is in agreement with our assumption concerning the 
in-plane coherence length £0 which is supposed to remain 
smaller than the size parameters of the SC grains [17]. 
A point of interest which should be addressed is the 
striking feature of the upper critical fields of the 
(TMTSF) 2 C104 measured along the a and the b direc- 
tions which are found to be greatly enhanced compared 
to the Pauli limit with an unusual upward curvature [21]. 



One may expect that the presence of SC and SDW do- 
mains in the (TMTSF)2C104 may act on the value of the 
critical fields as found in (TMTSF) 2 PFg characterized by 
the presence of SC slabs sandwiched by SDW phases [22] . 
A simple model based on this slab structure, which has 
been first proposed by Vuletic et at [3], gives rise to an 
enhancement and an upward curvature of the upper criti- 
cal fields of (TMTSF) 2 PF 6 [22]. However, we do not claim 
that the phase segregation picture proposed in the present 
work may explain the exotic behavior of the critical fields. 

Conclusion. — We have investigated, in this letter, 
the effect of the anion ordering on the SC transition tem- 
perature T c in the (TMTSF) 2 C10 4 salt taking into ac- 
count the superconducting fluctuations. The main result 
of this work is the presence of two regimes in the behavior 
of T c on quenching: a first regime with a rapid decrease 
for small cooling rates where SDW domains act as defects 
and a second regime with a smooth decrease ending by 
a saturation before the collapse of the SC phase at large 
cooling rates for which the SDW domains take over the SC 
ones. The obtained results are in good agreement with ex- 
perimental data. 

For a more accurate model, one may consider a random 
distribution of the SDW and SC islands. However, the 
general behavior obtained in the present work will be kept 
unchanged. 

Moreover, a complete description of the effect of the cool- 
ing rate on the SC phase should include the case where 
the size of the SC domains is smaller than the coherence 
length £0 as in granular superconductors. This peculiar 
issue goes beyond the scope of the present work. 
However, it is worth to stress that the phase segregation 
scenario, proposed in this letter, should be borne in mind 
to account for the striking features of the phase diagram 
of (TMTSF) 2 C10 4 . 
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